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Abstract 

A reanalysis of the radiative pion decay together with the calculation of the radiative 
corrections within chiral perturbation theory (CHPT) is performed. The amplitude of 
this decay contains an inner Bremsstrahlung contribution and a structure-dependent 
part that are both accessible in experiments. In order to obtain a reliable estimate of 
the hadronic contributions we combine the CHPT result with a large- A'^c' expansion and 
experimental data on other decays, which makes it possible to determine the occurring 
coupling constants. 



1 Introduction 



Rare decays are a useful source of information on particle interactions. Searches for new physics 
effects can take place at the high-energy or the high-precision frontier. At low energies new 
heavy particles can appear in the quantum loops. The advantage of high-precision physics is 
that one does not need to know the particle content of possible new physics in order to detect 
discrepancies between experimental results and theoretical predictions. One works with known 
external particles. On the other hand, it is of course not possible to detect new particles directly 
at low energies. 

The radiative pion decay tt"^ — > e^v^'y is interesting as it is not dominated by inner 
Bremsstrahlung and therefore sensitive to the so-called "structure-dependent" contributions 
that are generated by QCD effects. These are described with the help of two form factors, 
the vector and the axial-vector form factor. Via the conserved-vector-current (CVC) hypoth- 
esis the vector form factor can be related to the decays 7r° — >• 77 and 7r° — >• 7e+e~, where 
already precise data exist [1] . Therefore it is possible to measure the axial- vector form factor 
in experiments on the radiative pion decay directly. If one extracts both form factors experi- 
mentally with high precision, deviations from the CVC hypothesis can by investigated. Isospin 
breaking effects are of interest as they are not completely understood in the case of two-pion 
electroproduction and the corresponding r decay. There has been a discussion on a possible 
tensor interaction that could be detected in experiments on the radiative pion decay [2,3]. The 
induced tensor form factor due to radiative corrections is expected to be very small, but an 
explicit calculation seems to be useful. 

The typical energy scale of pion decays lies far below the region were perturbative standard 
model calculations are possible. At low energies chiral perturbation theory (CHPT) [4-7] is 
used as the effective field theory of the standard model in this energy region. All high-energy 
effects are included in the couphng constants of the effective Lagrangian. If it is not possible 
to determine all these coupling constants by use of experimental data from other processes, 
\axge-Nc QCD [8] is used to estimate them. 

Theoretical calculations of the structure-dependent contributions to the radiative pion de- 
cay have been presented in [9] and [10] at two-loop order. The lowest-order radiative corrections 
that are relevant for the inner Bremsstrahlung part are given in [11]. We complete the analysis 
by calculating the radiative corrections to the structure-dependent part within the framework 
of CHPT to lowest order in the large- A'^c expansion. 

On the experimental side an investigation has been performed at PSI [12], where the form 
factors have been measured with errors of only a few percent. Older experimental data [13-15] 
serve as an additional check. Future experiments on the radiative kaon decay will be helpful 
as the ratio of the decay widths of the pion and the kaon mode can be predicted theoretically 
with higher precision than the decay widths individually. 

The paper is organized as follows. In Sec. 2 we explain the basic facts of CHPT. The 
kinematics and the structure of the amplitude and the decay width of the radiative pion decay 
are presented in Sec. 3. The strong interaction contributions are explained in Sec. 4. Values of 
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the NNLO coupling constants are given. In Sec. 5 we report how the structure of the amphtude 
is modified by radiative corrections. Apart from that, the treatment of soft photon radiation 
and the apphcation of the Low theorem is explained and the large- A^^c form factors used to 
calculate the radiative corrections to the structure-dependent contributions are introduced. 
In Sec. 6 our results for the form factors, the radiative corrections and the decay width are 
presented. Our conclusions are summarized in Sec. 7. 



2 Low-energy expansion 

The asymptotic states in CHPT are not quarks and gluons but the members of the lightest octet 
of pseudoscalar meson^, the photon and the light leptons. CHPT is a low-energy expansion in 
the external momenta and masses that should be small compared to the natural scale of chiral 
symmetry breaking, which is expected to have a value of about 1.2 GeV. The order n of this 
expansion is indicated by p". A squared momentum of a pseudoscalar meson is of 0{p'^). The 
lowest-order effective Lagrangian is of the form 
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The pseudoscalar mesons are collected in a matrix u: 
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The symbol ( ) stands for the trace in three-dimensional flavour space. In order to introduce 
masses the external fleld x is set equal to an expression proportional to the quark mass matrix 
from now on: 

( rriu \ 
X = 2Bo\ . (2.4) 

V rris ) 

By adding terms determined via gauge symmetry to the external flelds and of the purely 
mesonic case the coupling of the photon and the leptons i, ui to the pseudoscalar mesons 



"'^ Other hadrons like the baryons can also be incorporated. 
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is fixed. 



r, = r^-eQ'^A^. (2.5) 

As the electroweak interactions break chiral symmetry, the spurion matrices Ql"r, Ql "^^^ be 
equated with the following expressions: 



Q 



L,R 



/ 2/3 


V 





-1/3 





Ql 



( Vud 


V 




(2.6) 



Every term in the Lagrangian (except kinetic terms) is multiplied with a coupling constant. 
The constant F in Eq. fl2.ip is identified with the pion decay constant in the chiral limit 
without electroweak interactions. In the same limit the constant Bq can be related to the 
quark condensate. Z dominates the pion electromagnetic mass difference. We use the SU(3) 
formalism in this work because information from processes involving strange quarks is needed 
to determine some of the coupling constants. 

The lowest-order Lagrangian is not enough to make connection with experiment. Higher 
orders have to be included. At every order the Lagrangian contains all terms that respect 
the symmetries. In [6] the SU(3) Lagrangian to O(p^) was presented considering also the 
Wess-Zumino-Witten functional [16]. Here and in the following only the terms relevant for our 
calculation are shown: 



\j+fivj+ J—fxuJ— I 
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I E^E^S^/^ + I + IT ^tiU^rM + Q{[va, vp])) + ... (2.7) 



where 



/r 
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Ud.UK 



(2.8) 



To 0(/) one has [18-20]: 

+C,2 {U,ufn + ^ C,s iU,, {x+, u'-u^}) + t C64 if+.uu'-un 

+Csi iUJn + Cs2 {U,u ir, X-]) + Cs7 {V,f-,.VT) 

+z ClTe^^"^ {x- \U,u. f-ap]) + CZe^"''^ MV,U,., + . . . (2.9) 
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with 



111 
= -[u\df,u] - -iu^r^^u - -iul^uK (2-10) 

The Lagrangian of 0{e'^p^) can be found in [21,22]. We will not present an 0{e^p'^) Lagrangian 
because at this order we use the expression that is of lowest order in the large- Xc expansion. 

As CHPT is a quantum field theory loops have to be taken into account. The primitive 
degree of divergence of a loop amplitude [17] is equivalent to the chiral dimension. A one-loop 
Feynman graph including only lowest-order vertices is of 0{p^) in the purely mesonic case and 
of 0{e^p^) if there is one internal photon propagator. The counterterms used to compensate 
the ultraviolet divergences of the loop integrals have to be of the same order in the external 
momenta as the loops. By renormalizing the appropriate coupling constants (e.g. Lg, C12, ...) 
that appear in the Lagrangian an UV finite amplitude is achieved. 

The scale dependent^] finite parts of all coupling constants are determined experimentally 
or estimated by performing resonance exchange calculations. In the large- iV^; limit the values 
of the coupling constants are given by exchange of infinitely narrow resonances. It turns out 
that at 0{p'^) the values one gets in this approximation agree quite well with the experimental 
values at a renormalization scale equal to the mass of the p particle [23]. This agreement is 
obtained by using only the lowest-lying vector, axial-vector, scalar and pseudoscalar octets. 

In [23] also the constant Z of Eq. fl2.ip has been determined. In this case and whenever 
one wants to calculate a coupling constant of an 0{e^p'^) Lagrangian resonance propagators 
appear in the loop and the correct momentum dependence of the involved form factors also 
for high energies [24] is needed. We summarize how this can be achieved in the case of the 
electromagnetic pion form factor F^. From Eqs. (12. ip and (12. 7p one get^: 

F,(g2) = 1 + ^ + 0{q') . (2.11) 

At leading order in the 1/Nc expansion including the lowest-lying vector resonance with mass 
My we have 

We will identify My with the mass of the p meson. The chiral loops are of higher order and 
introduce the width of the p [25]. Imposing that the form factor should vanish at infinite 
momentum transfer due to the Brodsky-Lepage behavior [26], the constant ky becomes equal 
to and 

Fe{q') = j^^,+0{l/Nc). (2.13) 



^The whole amplitude does not depend on the renormalization scale. 
^The renormalized coupling constants are labeled with an r. 
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Therefore, one concludes 



i«M^.) = 5j^. (2.14) 

The resonance Lagrangian that leads to a pion form factor with the correct low- and high- 
energy behavior to the order indicated in Eqs. fl2.1ip and (12.131) is of the form [24] 

(^res = W,,V.V^^ - ^ V,^Vn + ^ ( V/D + ^ {V,.U^U^) (2.15) 

with the high-energy condition FyGy = F^. The vector mesons are described by antisymmetric 
tensor fields V^u = ^ J2i=i KV^^. To 0{p'^) this formalism is equivalent to the more familiar 
notation with vector fields if one introduces explicit local terms [24]. 

3 General structure of amplitude and decay width 

The amplitude of 7r+(p) e'^{pe)i^{pu)l{k) has the following structure [27]: 

Mo = -zeGpV:,e;{F^L'^ - H'^'^Q (3.1) 

with 

= me^(p.)(l + 75)(^- 7 HPe), 

P ■ K ipe ■ K 

H^" = -^=^{Fv{pl)e^''^^k.pp-FA{pl){k-pg^''-p^k'')), 

F = U{p,h^{l-^^)v{pe), Pu,=Pe+Pu (3.2) 

where is the physical pion decay constant. One distinguishes between the inner 
Bremsstrahlung (IB) contribution and the structure-dependent (SD) part. The first is given by 
the term with and corresponds to the radiation of a pointlike pion and positron. The latter 
contains the two structure functions Fyipl^) and Fa(p^) including the hadronic contributions. 

In the process vr ewy the IB part is helicity suppressed, allowing the detection of the 
structure-dependent terms. The IB contribution diverges if the photon energy goes to zero. 
This divergence is canceled in the total rate by loop corrections to the decay ir ^ eu implying 
virtual photons. In experiments usually an energy cut is applied. Only photons above a certain 
energy are detected. 

Whereas the IB part is completely determined by the Low theorem [28] the structure- 
dependent part reflects the influence of QCD on this decay. The form factors Fy and Fa to 
(9(p^) in the chiral expansion are given by 

^ r?v+ 1 ^ ^ , ^ 

F^ 4v^7r2 ' ^ ^ 
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They include a mass 7717^+ that is of no physical meaning and drops out in the amplitude 
(see Eq. (13. 2p ). At higher orders and by including radiative corrections, the form factors get a 
momentum dependence. 

The importance of the different contributions can be seen from the differential rate (here 
normalized to the non-radiative mode): 



(l^r,.) = /i?(x,y)+f ^^^^-^ ) x[il + ^fSD^ix,y) + il-^fSD-ix,y)] 

(1 + 7) 5+ {x, y) + il- 7) Sr^ (a;, y)] (3.5) 



\V2F^ 

with 



^ = FA/Fy. (3.6) 

IB, SD^ and Stt are functions of the two kinematic variables x = 2p-k/m'l and y = 2p-pe/'m1. 
For iTie/mT^ = one has: 



IB{x,y) 



{l-ym + {l-xY) 



x'^{x + y — 1) 
SD+ix,y) = (l-x)(x + y-l)^ 

SD-ix,y) = (l-x)il-yf. (3.7) 

In Eq. (13.51) the terms including SD^ dominate over those with because of the additional 
factor m^/ml. When x + y goes to 1 the function IB{x,y) diverges. SD^{x,y) reaches its 
maximum at x = 2/3, y = I and SD~{x,y) at a; = 2/3, y = 1/3 (i.e. x + y = 1). One can 
define an angle between the positron and photon momenta: 

. 2 ^£7 X -\- y — 1 /o o\ 

sm — - = . (3.8) 

2 xy ^ ^ 

For Oe-y = the function IB{x,y) goes to infinity {SD^{x,y) has its maximum), whereas 
SD^{x,y) has its maximum for Oe-y = tt . Therefore an experiment performed in the region 
near O^^y = vr is sensitive to (1 + 7)^. It is difficult to distinguish experimentally between the 
terms proportional to IB and SD~ . 

In the standard model weak transitions are described by V— A interactions. New physics 
could lead to tensor interactions of the form 

T = ^^^^e;KFTu{p.)a'^^il + ^')v{p,) . (3.9) 
Radiative corrections generate an induced tensorial form factor as described in Sec. [5l 
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4 Contributions due to the strong interaction 

The form factors Fy and Fa have been calculated up to 0{j)^) for the chiral group SU(2) in [9] 
and for SU(3) in [10]. The momentum dependence of the form factors starts at 0{p^). We will 
use the SU(3) result which is in the isospin limit of the following form: 



AV27r^F^ 



3 ""-^^ + 3 + 327r2F2 



10 2 



with 



and 



G{z, rn? 



m 



1 - 



iln.^ 



4 + 



(4.1) 



(4.2) 



Fa{pI 



4^2 



(-2L^ + L^2Kln (^^^ 



<lng)+2.^1n(g; 



^ - 4v^ km\{QCl, - 2Q2 + ^64 + 2C8^i) 



6(27r)8 V 2m2 
+2m^(6C*{2 + 6Q3 — 2Cg^ — 2Cg2 + 2Cg3 + Cl^ + C^^ — C^g + 2C, 



80 



-2C, 



81 



2C*82 + C) 



87 



i(p^-m2)(2C,^3-4Q, + C8g 



(4.3) 



All coupling constants are taken at a scale equal to the p mass. For the two-loop integral I2 in 
Eq. (14. 3 p the numerical approximation given in [10] is used: 



Uz) = 44.5 z- 10304.2 



(4.4) 



The renormalized low-energy constants C^^, CI have to be determined by use of large- A'^c QCD 
or experimental data. Values for and can be obtained via the conserved vector current 
hypothesifl with the help of experimental data [1] on the decays 7r° 77 and 7r° — >■ 7 e^e~. 



exp 



\ 



2r(7rO 77) 



0.0262 ±0.0005. 



(4.5) 



The slope parameter of Fy is given by 



i!''P = 0.032 ± 0.004 . 



(4.6) 



^The relation that leads to Eqs. 
in case of the decay 7r° 76+6". 



and (|4.6p is reproduced within CHPT if one neglects the kaon loops 
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CUM,) 


Value [10-5] 


Source 


^12 


-0.6 ±0.3 


scalar resonance exchange 


/~ir 
^IZ 


0±0.2 


resonance exchange 


r^r 
'-'61 


1.0 ±0.3 


r decays, < W > correlator 


'-'62 


0±0.2 


resonance exchange 


^'-'63 '-"65 


1.8 ±0.7 


Kq charge radius 


'-'64 


0±0.2 


resonance exchange 


'-'78 


10.0 ±3.0 


resonance exchange 


/~ir 
'-"80 


1.8 ±0.4 


«!, Ki differences 


/^r 
'-'81 


0±0.2 


resonance exchange 


'-'82 


-3.5 ± 1.0 


resonance exchange 


/~ir 
'-^87 


3.6 ±1.0 


resonance exchange 


'-'88 


-3.5 ± 1.0 


resonance exchange 



Table 1: Values of the coupling constants appearing in Eq. (14. 3 p and the source of information 
used to fix them. 

One gets the following values for the low-energy constants in Eq. (14. ip 

"(Mp) = (0.1 ± 1.2) X lO-^MeV-2 , 
C^'{Mp) = (5.4±0.8) X lO-^MeV-2. (4.7) 

In [29] the constant C12 has been fixed by taking into account the exchange of scalar resonances. 
The constant Cei can be determined with the help of experiments on r-decays by considering 
the correlator of two vector-currents and using finite-energy sum rule techniques [30,31]. The 
combination 2Cq3 — Cq^ also appears in the expression for the electromagnetic Kq charge 
radius [32,33] that has been measured [34,35]. In [36] a large- iV^; expression for the correlator 
of vector, axial- vector and pseudoscalar currents with the correct high-energy behavior fixed by 
the operator product expansion is used to determine amongst others the low-energy constants 
C785 C*82, C*87 and Css- The contribution to 6*82 with three propagating resonances is not 
constrained by the high-energy behavior and will be neglected. The constant Cgo is fixed with 
the help of mass and decay constant differences of the ai and Ki particles following an idea 
presented in [31] for vector mesons (see App. [C]). The situation in the case of axial- vector 
mesons is more complicated as the states with the quantum numbers J^*" = l"*""*" and 1^^ 
mix. The other constants C13, Cq2, Cq4^, and Csi are set to zero as resonance exchange does 
not contribute in this case [37]. 

In Table [1] the values of the constants Cj at the p mass and the information needed for their 
determination are shown. 

The best procedure to get precise values for -Fy(O) and for the slope of this form factor is to 
take the values obtained in the isospin limit via the conserved vector current hypothesis (see 
Eqs. (14. 5 p and (14. 6p ) and to add/subtract the theoretical predictions for the isospin breaking 
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(ISB) contributions. The latter are proportional to m„ — m^, and 



In [38] the ISB contribution for 7r° ^ 77 has been calculated with the result 

Fv'^J'W = 0.00066 ± 0.0001 (2.5 % of F^ZViO)) . (4.9) 

What has to be considered concerning p^'^j^'^^'^^ ^re the radiative corrections discussed in 

ISB 

Sec. |5] and the following contribution proportional to an additional constant 

+ ^e+.7 ^ ^2^2 - ^q) 

From experimental data [39] on £+z/7 one gets| Cf[^(Mp) = (0.68 ±0.21) x lO^^MeV"^ 

which leads to f/^+^^^'^t = 0.00025 ± 0.00009 (0.9 % of F^ZViO)). Numerical results for the 

form factors can be found in Sec. El 



5 Radiative corrections 

The amplitude including radiative corrections contains additional terms compared to Eq. (13. 2p 
and is of the form 

M = -zGFeV:,e;{F^L'^FjBix,y)-H>^''Q+Tix,y), 
= -^—(e'^^-^V^^(x,y)-FAix,y)ik-pg^^-p^kn 

-FA{x,y){k-pig^^ -p^k^)), (5.1) 

where Vap^x.y) has a tensor structure more complicated then Fv{q'^)kaPf3 and one can dis- 
tinguish between two different axial-vector form factors. As mentioned above the radiative 
corrections generate in addition an induced tensorial form factor T{x,y), that is very small, 
i.e. 0.5 - 1.5 % of the IB part of the differential branching ratio depending on x and y. To 
O(e^p^) there is no contribution to Va/3{x,y) and the contributions to FA{x,y) and FA{x,y) 
turn out to be proportional to m^^/rn^ and can be neglected. Therefore the lowest-order radia- 
tive corrections (including also the induced tensorial form factor) can be regarded as corrections 
to the IB part. 

The squared amplitude |Mop receives radiative corrections due to virtual loop photons Ay 
and additional soft real photons A5: 

^ |M|' = ^ |Mo|' (1 + ^ (Ay + A5)) . (5.2) 

spin spin 

^We have assumed that the 0{p^) contribution is smaller then the 0(p'*) part. 
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Figure 1: Feynman diagrams of the virtual radiative corrections. The dot is a lowest-order 
vertex, the diamond is an vertex and the square can be an 0{p^) or an 0{p^) vertex. 

Diagrams due to wave-function renormalization are not shown. 
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Vertex 


expression 


7r+(t;)7r-(t)7*(g) 
7r+(i;)7r-(t)7(gi)7*(g2) 

vr+(t;)7(gi)7*fe)^+H 


^ (g^t;.e - v.qq.e) LI 

/I 2 

-^((ei.e2 ^2 - ei.g2 q2-e2)Ll + 2(gi.g2ei.e2 - gi.e2 ei.g2)(-^^9 + -^^lo)) 
-^{{(^w-f^q^ - e^.qq.e)Ll + {w.qe^.e - w.ee^.q){Ll + L{q)) 
^(g2-e«, ei-e2 - g2-ei e2-e«, + gi-e«, €1.62 - gi.e2 ei.e^)(L;; + L{q) 



Table 2: Even- intrinsic-parity vertices that are needed to calculate the loop amplitude of 
0{eY) 



Ac 



1 

'2^ 



P 



Pe 



(Pki 

2u;i \{kip) {kipe 



(5.3) 



aJi<A_B 



where tui is equal to \lk'l + X^ with a photon mass A introduced to deal with the infrared 
divergences and AE is the maximally allowed energy of the soft photon. The infrared divergent 
terms given by 



A5 

Ay 



inr^) + 




(5.4) 



cancel each other. The Feynman diagrams with a virtual loop photon are shown in Fig. [H 

The constant part of FiB{x,y) is obtained from the amplitude of tt'^ — > e~^u via the Low 
theorem: 



FjBix,y) = 1 + e' 
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-Ki + -K2 + — + —Kg + 2Ki2 
3 3 9 9 



'X-y — 2^2 



+2X3--X6 + Aioop(x,l/)). 



(5.5) 



The Ki and Xj are coupling constants of the 0{e^p'^) Lagrangian. By use of experimental data 
on the decay 7r+ fi^u no unknown low-energy constants remain to (9(e^p^). 

As the form factors Fy and Fa have been determined up to 0{p^) it makes sense to calculate 
also the radiative corrections of 0{e^p'^). In Table [2] we show the even-intrinsic-parity 0{p'^) 
vertices that are needed. The axial field is always replaced by 2v^ e+7^z/e^. The 

odd- intrinsic-parity vertices can be derived from Eq. fl2.7p . We will work at lowest order in 
the large- iVp expansion, therefore we will not include purely mesonic loops in our calculation. 
One also has to consider counterterms of 0{e^p'^). But as the coupling constants appearing in 
these counterterms are unknown we use large- A^^c form factors that include the 0{p^) vertices 
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and also produce the counterterm contributions as indicated in Eq. (12.121) in the case of the 
pion form factor. We will restrict ourselves to propagating p particles, for the ai a momentum 
independent contracted propagator will be used. By doing this one misses the contributions to 
the (9(e^p^) coupling constants coming from ai exchange which enlarges the error. We also do 
not consider resonance exchange in the odd-intrinsic-parity sector. 

It turns out that one just has to make the following replacements in the vertices in Table 
[2] in order to get the corresponding form factors obtained by use of the resonance Lagrangian 
Eq. (EISD. 

FyGy . 

Here qi is the momentum of the virtual photon. 

The sum of the loop graphs including the large- A^^c form factors is UV finite except the graph 
that contains the third vertex of Table [2] with an external photon and no loop photon. This 
divergence is due to the fact that the Oi propagator has been contracted. With a propagating 
ai in the loop also this graph is finite. So we make the following replacement in order to get a 
finite result: 



(47r)2 \pfl ) (47r)2 \d-A 2' ' ' y j ' jj ^^^y y^2^ 
Related to wave-function renormalization in the 0{e^p'^) amplitude, there is a term of the form 

(1 - ^ e\X^ - 47^12)) X 4v^m^+/F^ (L^ + L\^) (5.8) 

where F-,^ is the physical decay constant that includes also the 0{e^) contributions. As shown 
in (15.81) . we get a short-distance factor of the correct form [40] but it is unrenormalized. Our 
resonance calculation of the rest of the 0{e^p'^) amplitude that is finite by itself and valid to 
lowest order in the large- A^^c expansion does not generate the counterterm to renormalize the 
short-distance factor. One has to perform a two-step matching procedure of CHPT to Fermi 
theory and to the standard model [40] . We will use the following result [40] for the renormalized 
short-distance factor: 

Sew = l-le'iX;iM'^)-AK;,iM'^)) 

Putting together all contributions (see also App. [X]) and using the physical electron mass, the 
decay width with radiative corrections up to (9(e^p^) is of the form 

- Gl\V^,\'aSEw{'^^^^^^IBix,y)(l + -Ajsix,y)) (5.10) 



ax ay ^ 87r^ ^ tt 

;i + 7)' SD+ {x, y) + {l- 7)' SD- (x, y)]{l + -Asd{x, y)) } . 
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Figure 2: The relative size of the radiative corrections to the 0(e^p^) contribution (left figure) 
and to the 0(e^p'*) part (right figure) for y = 0.3, y = 0.5, y = 0.7 and y = 0.9. 

6 Results 

Without radiative corrections we predict the following form factors of Eq. fl3.2p : 

F^^^'^^"^ = 0.0262 ±0.0005 + (8.72 ± 1.09) X 10~^pym^ + O(pt), 

p^+^e+u^ = 0.0106 ± 0.0036 + (2.03 ± 0.65) x IQ-^pl/ml + 0{pi) (6.1) 

where is equal to m1{l — x). Except for the small isospin breaking contributions Fy is 
determined by data on the decays vr^ 77 and vr'' — 7e^e~. The 0{p^) contribution to 
Fa{0) is about 15 % and very sensitive to the values of the Lj. We have used the set of the Li 
given in App. [HI With the older set of values quoted in [10] the C(p^) contribution would be 
bigger. The relatively large error of -Fa(O) is due to the fact that the following sum of coupling 
constants is not known precisely: 

Li; + L^o = (1.39 ± 0.28) X 10"^ (6.2) 

In contrast to [10] we have quoted values for all of the appearing coupling constants and 
updated values are used. This is the reason for the difference of a few percent between the 
theoretical results presented in [10] and in this paper. 

In Fig. [2] we show the size of the radiative corrections ^ Ajb (left figure) and ^ ^sd (right 
figure) in percent depending on the kinematic variable x for four different choices of y. The 
first are negative over the whole phase space and smaller then 5 %. The latter are between -4 % 
and +4 %. The maximally allowed energy of the soft photon AE is set equal to 30 MeV. Up to 
a small difference that could be due to a misprint we agree with the result for Aib{x, y) in [11] 
under the assumption that in [11] the electron mass without radiative corrections and not the 
physical electron mass is used. The expression 1 + ^ AiB{x-,y) is used in [11] as an overall 
factor for the complete decay width and not only for the inner Bremsstrahlung part. This is 
only correct in the leading logarithmical approximation. But as all the dependence on mg/m^ 
cancels in the total decay width in accordance with the Kinoshita-Lee-Nauenberg theorem [41] 
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Tpmin 






Rthe 




(MeV) 


(MeV) 




(xlO"^) 


(xlO~^) 


50 


50 




2.58(8) 


2.655(58) 


10 


50 


40° 


14.77(40) 


14.59(26) 


50 


10 


40° 


38.89(90) 


37.95(60) 



Table 3: Theoretical (-Rthe) and measured (-Rexp) branching ratios for the three indicated phase 
space regions. 

one needs the CHPT expression to estimate the magnitude of the radiative corrections to the 
structure-dependent part of the total decay width. 

In Table O we compare the results for the branching ratios including all contributions with 
data [12] for experimental cuts indicated by E^'"', and 6'™". Our results agree with 

experimental data [12]. 

The fact that the theoretical branching ratio is very sensitive to Lg + L^g allows a rather 
precise determination of this sum of coupling constants if one uses the experimental result for 
the cuts = 50 MeV and E";''' = 50 MeV: 

{L,{Mp} + L,,[M^}} - I ^ Q ^g^ ^ ^Q_3 ^^^6) _ (6.3) 

This result is in good agreement with the existing theoretical prediction in Eq. (16.21) . 

7 Conclusions 

The radiative pion decay that includes an inner Bremsstrahlung part and a structure-dependent 
contribution has been reanalyzed. We have calculated the radiative corrections to the structure- 
dependent part to lowest order in the large- A^^c expansion within CHPT for the first time. 
Explicit values for all of the occurring 0{p^) coupling constants have been given by using and 
extending existing results [30-32,36]. 

It turns out that the 0{p^) contribution is about 15 %. The radiative corrections are 
a few percent varying over the phase space. The branching ratio agrees within the errors 
with experimental data [12]. There is no need to introduce a tensor interaction to explain 
the measured differential decay width obtained by use of the new data set [12]. The CVC 
hypothesis that relates the vector form factor of the radiative pion decay to the decay vr*^ 77 
seems to be a good approximation. 

The biggest theoretical error comes from the fact that the quite small sum of the coupling 
constants Lg and is not known with high precision. As a consequence a possible new physics 
contribution that affects the axial- vector form factor Fa is difficult to detect. Experimental 
data [12] allow a precise determination of Lg -|- L^q which also appears in the radiative kaon 
decays — »• e^Ugj and — > ^'^v^'^ and in Compton scattering 77?+ — > 771"'". 
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A Explicit form of the corrections Ajb and Asd 

The radiative corrections introduced in Eq. (15.101) are of the following form [11]: 

xix{y-l)-2y)\n\y) {{1 - y)x^ + 2{y - 2)x - Ay + A)\n{x)\n{y) 
iBK^^y) 4(x2-2x + 2)(y-l) (x2 - 2a; + 2) (y - 1) 

X (y^ + 1) ln(a; + y — 1) ln(y) x{x + y — 1) [y"^ + xy — 2y + x — \) ln^(x + y — 1) 



(a;2 - 2x + 2) (y - 1) 2 (x2 -2x + 2){x + y- 2)2 

\n{y) vr^ (-3(?/ - l)x'^ + 2{y'^ + ?,y ~h)x - 12{y - 1)) - 3 (x^ + 2) (y - 1) 



2 12 {x^ -2x + 2){y- 1) 

(jy - l)x' - 2{y - 2)x + 4(|/ - 1)) ln(l - x) ln(x) _ ^ / 2AE \ 
2(a;2-2x + 2)(y-l) ""XymJ 

+ In f In f tl!}jL \ + - + y - 1) + 2/ - 1) _ 3 / Z/'^^ 

\ymT, j \ ml ) {x'^ — 2x + 2) {x + y — 2) 4 \ . 

x{-yx + X + 2?/)Li2(l - x) {{y - l)x^ - 2{y - 2)x + A{y - 1)) Li2(x) 

^ 2 (a;2 - 2a; + 2) (y - 1) ^ 2 (x2 - 2a; + 2) (?/ - 1) 



xix + y~ l)Li2(l - y) ^ a;(x + y - l)Li2 (^) _ 3 / A£\ _ ^ 1 ^ ^^^^^ 



2(x2-2x + 2) 2(x2-2x + 2) 4 Vm2 y 2 

Aii^ is the maximal energy of the not detected additional soft photon and Ci , that is given 
by [22] 



1 , /M2\ Z 



+ In 



2 VM2 ) 4 



(A.2) 



has been defined in [42]. 



+ 



\ymT^J \ \ ml ) ) A \ ml J \ y 

l{x-l)iy- If Hy){32n^Ll + 32n^Ll, - if + ^^^^^^^ In i^) ((.(a; 
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+2y - 2) + 327r2(x^ + 2{y - l)x + 2{y - lf){Ll + L^o))« - A^)) 
+ ln ( -|- j (((x - l)(M2^(167r2(15x2 + 30(|/ - l)x + 13(|/ - lf)F'' 

+(25x^ + 50(y - l)x + 44(y - 1)^)M^^) - 16FV(15x^ + 30(y - l)x + 13(y - 1)^)M^^) 
+327r^L^(M„\(167r^(x - l)(15x^ + 30(7/ - l)x+n{y - lf)F'^ + (25x^ + (62y - 87)x^ 
+2( V - 55y + 46)x - 30(y - lf)Ml) - 16FV(x - l)(15x^ + 30(y - l)x 
+17(y - 1)^)M^2) + 327r^Lio(M„\(167r^(x - l)(15x^ + 30(?/ - l)x + 17(y - l)^)^^ 
+ (25x^ + (62?/ - 87)x^ + 2(V - 55?/ + 46)x - 30(?/ - lf)Ml) - 16FV(x - l){lbx'^ 

+30(y - l)x + n{y - lf)Ml)))/{l2MlMl) - ^^^U^^l - x)(l0247r^(x - 2){x^ 

+2{y - l)x + 2{y - lf)Ll^ + 6Att\32tt\x - 2){x^ + 2{y - l)x + 2{y - lf)Llo 

-x{x + 2y- 2))Ll + W2in\x - 2){x^ + 2{y - l)x + 2{y - lf)Lll - x{x^ + 2{y - l)x 

+2(y - 1)2) - 6A7r^x(x + 2y- 2)Lio) - \n(x) ((10247r^(x2 + Ax - A){x'^ + 2{y - l)x 

+2(y - lf)Ll'^ + 6A7r^(x(x^ + x - 2)(x + 2y - 2) + 327r2(x2 + Ax - A)(x'^ + 2{y - l)x 

+2{y - lf)Llo)Ll + 10247r^(x' + Ax - 4)(x' + 2{y - l)x + 2{y - l)^)Lll 

+ {x - 2)x{x'^ + 2{y - l)x + 2{y - if) + 6An^x{x^ + x-2){x + 2y- 2)L\q)') /{2x) 

+ l:ln{x + y - 1)( - 20A8n\x - l){3x^ + 6{y - l)x + A{y - if )Ll'^ - 32ti^{11x^ 
6 ^ 

+28yx^ - 39x^ + Ay^x - A2yx + 38x - m/ + 20y + 1287r^(a; - l){3x^ + 6{y - l)x 
+A{y - lf)L{Q - 10)L; - 20487r4(a; - l)(3a;2 + 6(y - l)x + A{y - iy)L[l 
-{x - l){5x^ + 10{y - l)x + 16(1/ - if) - 327r2(llx^ + {28y - 39)x^ + {Ay^ - A2y 
+38)x - 10{y - lf)Llo) + (l0247r^(x - l){6n'^{x'^ + x- 2){x'^ + 2{y - l)x + 2{y - if) 
+x{105x^ + 6(36?/ - 25)x2 + 2(31?/^ + 7y- 38)x + lAA{y - l)^))Ll'^M^^M^y^ 
+10247r^(x - l)(67r^(x2 + x-2){x'^ + 2{y - l)x + 2{y - 1)^) + a;(105a;^ + 6(36?/ - 25)x^ 
+2(31?/^ + 7y- 38)x + lAA{y - l)^))Lf^MlMy - Ut:^xL\^{AF^t:\x - l)x((3(y3 - 27y 
+18)^2 + 6(y^ - ?/3 - 27?/=^ + 45?/ - 18)x - (y - lf{92y^ + 81y - m^l, + {-3{by^ 
-27y + 18)x2 - 6(5?/^ - 5y^ - 27y^ + A5y - 18)x + {y - lf{20y^ + 81y - 54))M2) 
-y^{6Tr^{x + l){x + 2y - 2){x - 1)^ + x{96x^ + 4(56?/ - 71)x^ + (120?/^ - 439?/ + 283)a; 
-164y2 + 259y - 95))M^^M'^) - {x - l)x{8F^n^x{{3{y^ - 27y + 18)x^ + 6{y^ - y^ - 27?/^ 
+45?/ - 18)x +{y- 1)^(98?/^ - 81y + 54))M^^ + (-3(5?/^ - 27y + 18)x'^ - 6(5?/^ - 5y^ 
-27?/' + 45?/ - 18)x - (?/ - lf{50y^ - 81y + 5A))M^) - y'^iGn^x - l){x^ 
+2{y - l)x + 2{y - 1)') + x{87x^ + 6(25?/ - 24)x + 4(32?/' - 61y + 29)))M2^M^') 
+647r'L^(327r'(x - l)(67r'(x' + x - 2){x'^ + 2{y - l)x + 2{y - 1)') + x{105x^ + 6(36y 
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-25)x' + 2(3 V + 7y - 38)x + 144(y - iy))L\^Ml^My + x{y\ioT:\x + l)(x + 2y 
-2)(x - 1)^ + x(96x3 + 4(56y - 71)x^ + (120y^ - 439y + 283)a; - 164y^ + 259y 



95))M4 Af^ - 4FV^(x - l)x((3(y3 - 27y + l%)x' + 6(y^ - y 



-{y - l)^(92y3 + 81y - 54))Mf^ + (-3(5?/3 - 21y + 18)x^ 



-1 



,x + (y - l)'(20|/3 + 812/ - 54))M2))))/(36xVM^, M, 
-1)2 + (x + y - l)2(327r2(L^ + L\,) + 1)^) (1 - x) 



27y2 + 45y - 18)x 
6(5y^ - - 27y2 + 45?/ 



I Ml 



(A.3) 



B Numerical input 

In this appendix we collect the numerical values of coupling constants and masses used in this 
article that have not been already explained before. 

Masses [1] 



Mp = 775 MeV 



139.570 MeV 



= (m^o + m7r+)/2 



m^o = 134.977 MeV 



Chiral low-energy constants [40,43-45] 



F = 87.7 ±0.3 MeV 
LI = (0.43 ±0.12) X 10-3 
(6.49 ±0.20) X 10-3 
(0.7 ±0.3) X 10-3 
= (-4.2 ± 1.5) X 10 



LI 



-3 



F^ = 92.2 ± 0.3 MeV 
= (0.73 ±0.12) X 10-3 
-5.10 ±0.20) X 10-3 



Tr 



(11.6 ±3.5) X 10- 



^6 



-2.56 ±0.50 
(-2.53 ±0.37) X 10-3 
(-2.7 ±0.9) X 10-3 
(2.8 ±0.9) X 10-3 



Concerning Lg and L[q we have used the mean value of the following determinations: 



5.93 ± 0.43 X 10-3 [46] 
7.08 ±0.40 X 10-3 [23,36] 
-- 6.46 ±0.40 X 10-3 [47,48] , 



2M2 

_5 ]_ 

2V6 167r2 



(B.l) 



-5.13 ±0.19 X 10 



15 1 
8^6 167r2 



[49] 

-5.31 ±0.40 X 10-3 [23,36] 
= -4.85 ± 0.40 X 10-3 [47,48] . 



(B.2) 
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C Determination of CgQ 

The constant CgQ in Eq. (14.31) can be determined via resonance saturation by use of mass and 
decay constant differences of the axial- vector mesons ai and Ki. The relevant terms of the 
resonance Lagrangian [37] are 

Cr = -^{A,,An + \t^{x+A^uAn + >^'''''{SA,,An (C.3) 
+^{A,jn + ({5, A,^}fn - ^MliS') + c,{Su^u,) + c^(%+) 

where the antisymmetric tensor field A^^, contains the axial-vector mesons and S includes the 
scalar mesons. One obtains the following expression for Cgo [37]: 



(^r _ p2 ( CdCm ( ^AA _ 9 /o\SA ^ACm ^SAA -^j^m 

In analogy to the notation in [31] we define: 



(C.4) 



,AA I ^rn_ ^SAA 
\SA 

The physical /ri(1270) and i^'i(1400) states are a mixture of the J^*" = l^"*" and 1+^ states 
KiA and Kib'- 



eA=2 Ur + ^A^--) , (C.5) 



fTi = Ma,j^Xl\ (C.6) 



fs:i(1270) = KiAsm9 + KiB cos9, 

Ki{UOO) = K^A cose -KiB sine. (C.7) 
With a mixing angle e of (59 ± 3) degree [50] the mass of the Kia state is given by 

Mk,^ = 1308 ± 10 MeV . (C.8) 
There is the following relation [51,52] between the decay constants of the Kia and the i^'i(1270): 

Fk^^ = . ^ = 159 ± 20 MeV (C.9) 
sm b' — cos e 

with [51] 

1 m, — m,, ,^ 

6 = ^— ^^0.16. C.IO 

V2 rris + rriu 

One obtains values for and fAi via the relations 

Mk,^ - MP'^ = AeJBo{ms - (m„ + m,)/2) (C.ll) 
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and 

Fk,, - = ^-^i?o(m, - (m. + m,)/2) . (C.12) 

Together with Fa, = 165 ± 13 MeV [43], M^f = 1230 ± 40 MeV and m, = 25.90 m„ this leads 
to 

Cg^o = (1.8 ± 0.4) X 10-\ (C.13) 
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